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ABSTRACT 



Aims. We exploit deep observations of the GOODS-N field taken with PACS, the Photodetector Array Camera and Spectrometer, on 
board of Herschel, as part of the PACS Evolutionary Probe guaranteed time (PEP), to study the link between star formation and stellar 
mass in galaxies to z ~ 2. 

Methods. Starting from a stellar mass - selected sample of ~ 4500 galaxies with mag 45 ^ m < 23.0 (AB), we identify ~ 350 objects 
with a PACS detection at 100 or 160 pm and ~ 1500 with only Spitzer 24 pm counterpart. Stellar masses and total IR luminosities 
(L IR ) are estimated by fitting the Spectral Energy Distributions (SEDs). 

Results. Consistently with other Herschel results, we find that L IR based only on 24 pm data is overestimated by a median factor ~ 1.8 
at z ~ 2, whereas it is underestimated (with our approach) up to a factor ~ 1.6 at 0.5 < z < 1.0. We then exploit this calibration to 
correct L tR based on the MlPS/Spitzer fluxes. These results clearly show how Herschel is fundamental to constrain L /R , and hence the 
SFR, of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of 
the relations between the star formation rate (SFR), the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, 
in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR 
information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, 
being almost flat at z < 1 .0 and reaching a slope of a = -0.50+jj |g at z ~ 2, at odds with recent works based on radio-stacking analysis 
at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ~ 15 for massive M > 10 n M o galaxies from 
z = to z = 2, and seems to flatten at z > 1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, 
implying that they have formed their stars earlier and more rapidly than their low mass counterparts (downsizing). 

Key words. Galaxies: evolution - Galaxies: active - Galaxies: starburst - Cosmology: observations - Infrared: galaxies 



1. Introduction 

The link between galaxy stellar mass and star formation rate 
(SFR), and its cosmic evolution is crucial to shed light on 
the processes of galaxy formation. The specific SFR (SSFR = 
SFR/mass) plays an important role as it measures the star for- 
mation efficiency of a galaxy and the fraction of a galaxy mass 
can be converted into stars per unit time (see e.g. de Cunha et al. 
2010). Several studies at < z < 3 report similar findings: (1) 
the SSFR increases with redshift at all masses; (2) the SSFR of 
massive galaxies is lower at all z (e.g. Feulner et al. 2005, Erb et 
al. 2006, Perez-Gonzalez et al. 2008, Damen et al. 2009, Dunne 
et al. 2009). However, the scatter and the exact slopes of these 
relations are still not clear. In particular, the dependence of SSFR 
on mass is one of the most debated open questions. On the one 
hand, radio-stacking analysis of galaxies at z — 1 .5 -2.0 selected 
in the T^-band found a weak (or even absent) SSFR-mass correla- 
tion, with an indication of steepening at higher redshift (Pannella 
et al. 2009, Dunne et al. 2009). On the other hand, other results 
based on UV to mid-IR SFR tracers indicate a clear decrease of 
the SSFR with increasing mass (e.g. Feulner et al. 2005, Erb et 
al. 2006, Noeske et al. 2007b, Cowie & Barger 2008). Some of 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



these discrepancies could be due on one hand to the effects of 
dust extinction corrections for the SFR estimates based on UV- 
optical indicators, and on the other hand to the assumptions and 
large extrapolations on the infrared SED shape and luminosity 
adopted for the SFR estimates based on 24/im data. 

The advent of the Herschel Space Observatory (Pilbratt et 
al. (2010) finally allows us to robustly derive the total infrared 
(IR) luminosity (L/r) of galaxies, by directly sampling the peak 
of the thermal emission of dusty galaxies up to z ~ 3. In or- 
der to place new and stringent constraints on the SSFR evolu- 
tion, in this work we exploit the PACS Evolutionary Probe (PEP) 
guaranteed time data collected in the GOODS-North field with 
the PACS (Poglitsch et al. 2010) instrument. The PEP observa- 
tions are described in Berta et al. (2010, this Issue Appendix 
A). Ancillary data, both photometric and spectroscopic, in- 
cluding UV (GALEX), optical (HST), near-IR (FLAMINGOS, 
IRAC/Spitzer) and mid-IR (MlPS/Spitzer) data, have been col- 
lected to build up a reliable multiwavelength catalog and photo- 
metric redshifts. We adopt h - 0.7, CIa = 0.73 and Q m = 0.27. 

2. Sample selection 

Our aim is to investigate the evolutionary link between stellar 
mass and star formation avoiding strong biases and selection 
effects. Thus, the main galaxy sample was selected at 4.5/Km 
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Fig. 1. Top panel: Comparison of L/« derived from the UV-to- 
PACS/im data with the Polletta templates and the Chary & Elbaz 
(2001) libraries, by fitting one single point (100/im). Bottom 
panel: Comparison of L/« derived from the UV-to-24yum and 
the UV-to-PACS SED fitting procedures. Different colors refer 
to various redshift intervals (see legend in the figure). The thick 
cross symbols represent the median value in the corresponding 
luminosity bins. We also report with dotted lines the +l<x limits 
around the median. In both panels only sources from the 4.5 pm 
catalog that were detected with SNR> 3 at all 24, 100 and/or 
160 pm are reported. 



with IRAC (mag4 5 /jm < 23.0, AB) in order to ensure sensitiv- 
ity to stellar mass up to z ~ 2 - 3. At this limiting magni- 
tude, the IRAC sample is ~80% flux complete (Mancini et al. 
2009), it is not strongly affected by confusion (see Rodighiero 
et al. 2010) and it includes 4459 sources. On the IRAC po- 
sitions we fitted PSFs to the MIPS and PACS images (how- 
ever, for detections in the PACS maps we used only IRAC posi- 
tions with a 24um detection). Out of the 4459 IRAC sources, 
1887 (351) sources have MIPS (PACS) fluxes with signal-to- 
noise ratio (SNR) greater than 3. The typical 24pm, 100;um and 
160^m fluxes in this catalog reach the 3-<x limit, that is ~20^Jy 
(Magnelli et al. 2009), ~3mJy and ~5.7mJy, respectively. With 
this selection we miss only 2 PACS-detected objects with low 
signal-to-noise ratio (SNR~ 3). About 40%, ~ 52% and ~ 70% 
of the IRAC, MIPS and PACS subsamples have an optical spec- 
troscopic redshift (mainly from Barger et al. 2008), respectively. 



3. Stellar masses and IR luminosities 

Following the approach of Rodighiero et al. (2007), stellar 
masses were estimated by setting the redshift (photometric 
or spectroscopic) of each IRAC-selected object and using the 
Hyperz code (Bolzonella et al. 2000) applied to the photomet- 
ric SEDs in the optical-to-5.8^m range. For an easier compari- 
son with literature data, we used the stellar-population synthe- 
sis models of Bruzual & Chariot (2003, hereafter BC03) with a 
Salpeter IMF, exponentially declining r models for the SFR, and 
solar metallicity. We estimated the stellar mass completeness as 
a function of redshift as described in Mancini et al. (2009) (see 
online Figure lB~TT i. 



For the objects detected with PACS and MIPS, the IR lumi- 
nosities L/r were estimated with a best fitting procedure by com- 
paring the optical-to-IR SEDs with a library of template SEDs of 
local objects from Polletta et al. (2007) and adding a few modi- 
fied templates (see Gruppioni et al. 2010). We followed the same 
approach that we applied to 24pm GOODS-South and SWIRE- 
VVDS sources in Rodighiero et al. (2010) and in Gruppioni et 
al. (2010). The inclusion of the whole SED in the fitting pro- 
cedure allowed us to fully exploit the photometric information 
and determine the K-correction in the most reliable way. As de- 
scribed in Rodighiero et al. (2010), we forced the spectral fit 
to reproduce the far-IR Spitzer-Herschel datapoints by reduc- 
ing their photometric errors. To minimize the contribution from 
AGN emission to the L/«, we removed from our sample X-ray 
detected sources with Lx > 10 42 erg/s and those objects clas- 
sified by the SED fitting analysis as Type 1 quasars. A few ex- 
amples of our best-fits are presented in the online Appendix C, 
Figure ICTl 



4. The estimate of SFR 

For the galaxies detected with PACS or at 24pm, the instan- 
taneous SFR was estimated using the combination of IR and 
UV luminosity as in Papovich et al. (2007) and Santini et al. 
(2009): SFR IR+uv /M yr l = 1.8 x lO" 10 x L ho ,/L with L ho , = 
2.2 X Ljjy + L[r, where L/« was derived by integrating the best-fit 
SEDs in the [8-1000]//m rest-frame range. The rest-frame UV 
luminosity (which accounts for for the contribution for young 
unobscured stars), uncorrected for extinction, derived from the 
SED fitting with BC03 model, is L uv = 1.5 x L(2700A). The 
contribution from Luv is marginal, and its exclusion does not 
change the results presented here. In the rest of the paper, we 
will then discuss the main evolutionary trends linking mass and 
SFR by considering only A >24pm detections, and we will dis- 
cuss the contribution of the remaining far-IR undetected IRAC 
sources through a stacking analysis. 



5. Results 

5. 1 . Spitzer vs Herschel IR luminosities 

The reliability of our IR luminosities was assessed by compar- 
ing them with alternative approaches. For PACS detections, we 
compare in Figure Q] (top panel) our L/« derived from SED fit- 
ting to the UV-PACS range with Polletta templates, with those 
computed by using the Chary & Elbaz (2001, hereafter CE01) li- 
brary and fitting only the 100/im fluxes. This shows a very tight 
correspondence, confirming the value of far-IR information to 
constrain the bolometric energy budget of high-z galaxies, inde- 
pendently of the models assumed. For sources with 24pm detec- 
tion and no PACS, we have to rely on large extrapolations, that 
can be strongly influenced by the adopted models. As an inter- 
nal check, we then compare in Figure Q] (bottom panel) the IR 
luminosity derived from our UV-to-24/im and the UV-to-160/im 
SED fitting procedures. Different colors refer to various redshift 
intervals. As already discussed by Elbaz et al. 2010 and Nordon 
et al. 2010, at z > 1 .5 L« based on 24pm data are overestimated, 
by a median factor ~ 1 .8 in our analysis: this factor is larger for 
the CE01 models and when including PACS upper limits in the 
analysis (e.g. up to a factor ~ 4, see Nordon et al. 2010, showing 
that the observed trend is not induced by any PACS incomplete- 
ness). At z < 1.5, the 24pm based L/r are underestimated, up 
to a factor ~ 1.6 around z ~ 0.5 - 1.0. This result highlights 
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Fig. 2. Relation of the stellar mass as a function of SFR (top panels) and SSFR (bottom panels). Green points represent PACS 
detected sources for which the L/« has been computed including all data points. Black points are MIPS 24/im sources undetected by 
PACS: in this case the Ljr has been derived from our SED fitting procedure limited at 24/im. Open symbols are passive sources with 
color mag[24]-mag[3.6]<0.5. We also report with red lines the average of the relations and the ±l<x limits. The red dot-dashed lines 
correspond to literature observed relations between mass and SFR: at z ~ 1 the GOODS relation by Elbaz et al. (2007); at z ~ 2 
the relation for star-forming BzK by Daddi et al. (2007). Vertical dashed lines mark the mass completeness limit as a function of z, 
and the horizontal dotted lines indicate the inverse of the age of the universe at the mean redshift of each bin. The dashed blue lines 
in the bottom panels show the relations derived from our stacking analysis on PACS maps. Open circles mark sources classified 
as passive on a color basis (mag[24]-mag[3.6]<0.5), and generally populate a separate sequence (this bimodality has already been 
detected by Santini et al. 2009). 



the key role of Herschel and have strong implications on pre- 
vious Spitzer works on the SFR of z ~ 2 galaxies. In the fol- 
lowing analysis and figures we have corrected the 24//m based 
Lir, according to the new Herschel calibration that we have just 
discussed. 



5.2. The mass versus SFR relation of star-forming galaxies 

The existence of a strong correlation between galaxy SFR and 
mass at different redshifts (0 < z < 3) has been extensively 
discussed in the recent literature (Erb et al. 2006, Noeske et al. 
2007a, Elbaz et al. 2007, Daddi et al. 2007, Perez-Gonzalez et 
al. 2008, Pannella et al. 2009, among the others), but there is not 
yet agreement on the general validity and properties of this re- 
lation. The main results of our analysis are shown in Figure [2] 
where we combine the subsamples of the star-forming galaxies 
detected with PACS and/or MIPS. The upper panels indicate the 
existence of a (rather scattered) positive correlations between the 
SFR and stellar mass at all redshifts. The Spearman test indicates 
that the probability of not having correlation is 2e-23, 7e-23, 9e- 
41, le-28, 2e-14, 2e-17, 4e-08, and le-10 in the eight redshift 
intervals, respectively. The SFR-mass relation looks more robust 
at z < 1. However we should remind that in this case we are 
considering also sources below the mass completeness (vertical 
lines), and that we are sampling only the most star-forming pop- 
ulation, as MIPS or PACS undetected sources are not included 
in these plots. The comparison with Elbaz et al. (2007) at z ~ 1 
and Daddi et al. (2007) at z ~ 2 (see caption of Figure|2]for de- 
tails) shows that their observed slope of the SFR-mass relation is 
not inconsistent with our results, but even much flatter relations 
would be allowed by our data above the completeness limit. The 
lower panels of Figure [2] show the relation between SSFR and 
stellar mass for the same sources plotted in the upper panels. A 
negative trend of SSFR with mass is evident at all redshifts, al- 
though the scatter is quite large. The bulk of PACS and/or MIPS 
sources is located above the horizontal dotted line (the inverse of 
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Fig. 3. Relation of the stellar mass as a function of the SSFR for 
PACS detected sources in various redshift bins: the SFR for red 
points has been computed from PACS fluxes, while for black 
points it has been extrapolated from the SED fitting from the 
24yimi data. 



the age of the Universe), indicating that these systems are expe- 
riencing a major episode of star formation, forming stars more 
actively than in their recent past and building up a substantial 
fraction of their final stellar mass. 

The key role of Herschel for high-z galaxies is highlighted 
in Figure [3] showing that when PACS data are used, the scatter 
of the SSFR - mass relations decreases by factors of 1.0, 1.3, 
2.3 and 3.7 for increasing redshift bins, compared to the case 
of SSFRs estimated with MIPS data only. We have to warn that 
PACS detects only the brightest objects and we cannot than ver- 
ify that the scatter is intrinsically lower. However, the fact that 
at least at high luminosities, at z ~ 2, PACS produces a smaller 
scatter (because it provides a more accurate SFR), might suggest 
that a similar trend should happen also at low luminosities. 
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6. Stacking analysis on PACS images 

To test if the results in Figure 2 might have been affected by 
selection biases and to verify their reliability, we performed a 
stacking analysis (Bethermin et al. 2010) including all sources 
of the original IRAC [4.5]<23.0 sample. We splitted the sam- 
ple in bins of mass and redshift, and stacked on a residual 
160/im map (created by removing all PACS 160pm detections 
included in our catalog with SNR> 3) at the positions of all 
IRAC sources undetected by PACS (stacking at lOO/vm does 
not change our results). With this procedure, we derived a typ- 
ical flux for each mass-redshift bin. Using the formalism in- 
troduced by Magnelli et al. (2009), that accounts both for de- 
tections and no-detections, we then converted these fluxes into 
luminosities by adopting an average ^-correction for each red- 
shift (derived from the more common SED template adopted by 
our best-fitting procedure in that interval). An empirical relation 
between rest-frame vL v (l60pm) and Lm was used to first con- 
vert the stacked luminosities into bolometric luminosities, and 
then into SFR through a standard Kennicutt law (Kennicutt et 
al. 1998). In this analysis we have included only star-forming 
galaxies. To exclude passive sources we applied an empirical 
color selection of ( U - B) r f < 1.1, calibrated from our data, and, 
to recover massive dusty sources that might fall into the red se- 
quence, we included in the stacking analysis also sources with 
(U - B) rf > 1.1 andmag[24]-mag[3.6]>0.5. 

The results of the PACS stacking analysis are presented in 
Figure |4] (upper panel) and compared with literature data based 
on star-forming samples (with the exception of Dunne et al. 
2009 that includes all /f-band selected objects). The slope of our 
SSFR-mass relation becomes steeper with redshift, being quite 
flat at z < 1 with a - -0.25 *5 - H (all the slopes have been com- 
puted in a mass complete range). Intringuingly, the slopes are 
consistent with those observed for the 24, 100 and 160jum de- 
tections at 0.5 < z < 2 (blue, long-dashed lines in Figure 3). 
At z < 1, our results are in broad agreement with those based on 
radio-stacking that found almost flat relations up to z ~ 2 (Dunne 
et al. 2009, Pannella et a. 2009) , while at z > 1 our relation 
evolves toward stronger dependencies (a = -0.50^'jj?). We note 
that the rejection of passive and massive sources in the Dunne 
et al. sample would produce even flatter relations, increasing 
the discrepancy with our results at high-z. At low z, the study 
by Noeske et al. (2007) based on the UV-optical luminosities is 
steeper then our results, while at z ~ 2 the work by Erb et al. 
(2006) is more consistent with our observed SSFR-mass relation 
(however, within a very huge scatter). A recent work by Oliver 
et al. (2010), based on MIPS 70 and 160yum stacking, seems in- 
stead to be in general agreement with us up to z < 1.5, where 
they have enough statistics. 

Figure 4 (bottom panel) also shows that the mean SSFR of 
star- forming sources rises with redshift, up to a factor ~ 15 for 
the most massive galaxies (M > 1O U M ), implying that galaxies 
tend to form their stars more actively at higher redshifts. The 
mean SSFR seems also to flatten at z > 1.5 for (M > 1O 1O5 M ). 
Moreover, the most massive galaxies have the lowest SSFR at 
any redshifts. As shown in bottom Fig. 4, they have already so 
large stellar masses at z=0.7 to 2.5 that they would require steady 
SFR at the PACS observed level during the whole Hubble time 
at that redshift to form. 
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PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI/INAF (Italy), 
and CICYT/MCYT (Spain). 



References 

Barger, A. J., Cowie, L. L., & Wang, W.-H. 2008, ApJ, 689, 687 

Berta et al., 2010 this volume 

Bethermin, M., Dole, H., Beelen, A., & Aussel, H. 2010, arXiv:1001.0896 

Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 363, 476 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Cowie, L. L., & Barger, A. J. 2008, ApJ, 686, 72 

Daddi, E., et al. 2007, ApJ, 670, 156 

da Cunha, E., Eminian, C, Chariot, S., & Blaizot, J. 2010, MNRAS, 273 

Damen, M., Labbe, I., Franx, M., et al. 2009, ApJ, 690, 937 

Dunne, L., et al. 2009, MNRAS, 394, 3 

Elbaz, D., et al. 2007, A&A, 468, 33 

Elbaz et al., 2010 this volume 

Erb, D. K., Steidel, C. C, Shapley, A 

Feulner, G., Gabasch, A., Salvato, M. 

Gruppioni et al., 2010 this volume 

Kennicutt, R. C, Jr. 1998, ApJ, 498, 541 

Magnelli, B., Elbaz, D., Chary, R. R., et al. 2009, A&A, 496, 57 

Mancini, C, Matute, I., Cimatti et al. 2009, A&A, 500, 705 

Noeske, K. G., et al. 2007, ApJ, 660, L43 

Nordon et al., 2010 this volume 

Oliver, S., et al. 2010, arXiv: 1003.2446 

Pannella, M., et al. 2009, ApJ, 698, LI 16 

Papovich, C, Rudnick, G.; Le Floc'h, E. t al. 2007, ApJ, 668, 45 



E., etal. 2006, ApJ, 647, 128 
et al. 2005, ApJ, 633, L9 



G. Rodighiero et al.: The mass-SFR link in high-z galaxies 



Perez-Gonzalez, P. G., Trujillo, I., Barro, G, et al. 2008, ApJ, 687, 50 

Pilbratt et al., 2010 this volume 

Poglitsch et al., 2010 this volume 

Polletta, M., et al. 2007, ApJ, 663, 81 

Rodighiero, G, Cimatti, A., Franceschini, A. et al. M. 2007, A&A, 470, 21 

Rodighiero, G, et al. 2010 A&A preprint doi 

http://dx.doi.org/10.1051/0004-6361/200912058 
Santini, P., et al. 2009, A&A, 504, 751 



G. Rodighiero et al.: The mass-SFR link in high-z galaxies, Online Material p 1 



Appendix A: Authors affiliations 

1 Department of Astronomy, University of Padova, Vicolo 
dell'Osservatorio 3, 1-35122 Padova, Italy, e-mail: 
giulia.rodighiero@unipd.it 

2 Dipartimento di Astronomia, Universita di Bologna, via 
Ranzani 1, 1-40127 Bologna, Italy 

3 INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, 
1-40127, Bologna, Italy 

4 Max-Planck-Institut fiir extraterrestrische Physik, Postfach 
1 3 1 2, 8574 1 Garching, Germany 

5 European Southern Observatory, Karl-Schwarzschild-Str. 2, 
85748 Garching, Germany 

6 INAF-Osservatorio Astronomico di Roma, via di Frascati 
33, 00040 Monte Porzio Catone, Italy. 

7 ESAC, Villafranca del Castillo, ES-28691 Madrid, Spain 

8 Laboratoire AIM, CEA/DSM-CNRS-Universite Paris 
Diderot, IRFU/Service d'Astrophysique, Bat.709, CEA- 
Saclay, 91191 Gif-sur-Yvette Cedex, France. 

9 Instituto de Astrofisica de Canarias, 38205 La Laguna, 
Spain. 

10 Department of Astronomy, 610 Space Sciences Building, 
Cornell University, Ithaca, NY 14853, USA 

11 INAF-IFSI, Via Fosso del Cavaliere 100, 1-00133 Roma, 
Italy 

12 INAF-Osservatorio Astronomico di Padova, Vicolo 
dell'Osservatorio 2, 1-35122 Padova, Italy 

13 INAF-Osservatorio Astronomico di Trieste, via Tiepolo 11, 
34143 Trieste, Italy 
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Fig. B.l. Mass-completeness thresholds as a function of redshift 
for our IRAC 4.5/im-selected sample (mag4 5 A , m < 23.0, AB), 
derived from synthetic stellar population models as described 
in Mancini et al. (2009). Here we used the constant SFR tem- 
plates of Bruzual & Chariot (2003), with a Salpeter IMF, and 
different ages, and dust extinction parameters (Eb-v)- For each 
specific age, we considered three possible values of dust ex- 
tinction: Eb-v - 0.3 (solid lines), Eb-v - 0.5 (dashed lines), 
and Eb-v = 0.8 (dot-dashed lines). In our analysis we adopted 
the most conservative mass-completeness limit (dot-dashed ma- 
genta line), above which even the oldest (2 Gyr) and highly ex- 
tincted (Eb-v - 0.8) star-forming galaxy population would be 
entirely recovered. 
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Fig. C.l. The observed Spectral Energy Distributions (filled circles) of four PACS sources from our sample, representative of the 
various spectral classes considered in this work. We also show the best-fit spectra obtained with Hyperz based on the spectral library 
by Polletta et al. (2007) and the updated version of Gruppioni et al. (2010): the red solid lines represent the best-fit obtained with 
Herschel data, while red dotted lines are the fit limited to the 24/mi Spitzer data points. The blue lines show the fit to the optical-5.8 
yum bands with Bruzual & Chariot (2003) models . 



